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Abstract 

We verify the Edgeworth expansion of any order for the integrated ergodic Levy driven 
Ornstein-Uhlenbeck process, applying a Malliavin calculus with truncation over the Wiener- 



(N 



Poisson space. Due to the special structure of the model, the coefficients of the expansion can be 



< 

& 

■i— > 



given in a closed form. 
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1 Introduction 

Let (X, Y) = {(X t ,Y t )} t £R + be the bivariate model described by 

X t = X Q - A f X s ds + Z u 

t Jo 
Y t = f (7 + pX s )ds + pZ u 
Jo 



(1) 



, £». | where Z — (Zt)teR+ is a non-trivial Levy process independent of the initial variable Xq, and the 

iy-v . parameter (A, 7, /3, p) G (0, 00) x R x (K\{0}) x K satisfies that 

(N i P + M ? 0. (2) 

pi J The process X is the exponentially ergodic Levy driven Ornstein-Uhlenbeck (OU) process; we refer to 

(f-) . [4] and the references therein for fundamental facts concerning the OU process. The goal of this note 

is to provide conditions under which the Edgeworth expansion of the expectation £ , [/(T _1 / 2 iJ'r)] as 

T — > 00 is valid, where 

H T ■= Y T - E[Y T ] (3) 

^S ■ 

and / : R — > R is a measurable function of at most polynomial growth. The condition (2) will turn 

to be necessary for the Gaussian limit of C(T~ 1 / 2 Ht) to be non-degenerate: as a matter of fact, the 

necessity of (2) can be seen concisely by the expression 

T- l ' 2 H T = (/3 + pX)T- 1 ' 2 f (X t - E[X t ])dt + P T- 1 ' 2 {(X t - E[X t ]) - (X - E[X })} , 

so that, if (3 + p\ = and (X t - E[X t \) - (X - E[X ]) = O p (l) as T -> 00, then C{T- l / 2 H T ) tends 
in probability to (See Section 2.2). 

As is well known, distributional regularity of the underlying model is essential to the validity of 
the Edgeworth expansion. At first glance, the regularity of the joint distribution C(X, H), which will 
play an essential role in derivation of the expansion (see Section 3), does not seem enough since we 
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have only one-dimensional random input Z against the two-dimensional objective (X,H). In partic- 
ular, for pure-jump Z we have to take distributional regularity over the Poisson space into account, 
rendering the problem mathematically interesting in its own right. In this case, we will execute the 
Malliavin calculus under truncation, which enables us to successfully pick out a nice event on which 
the integration by parts formula can apply to ensure distributional regularity; more specifically, our 
truncation functional will be constructed through two diffusive jumps, so as to make the Malliavin 
covariance matrix associated with the flow of (X, H) non-degenerate (As will be mentioned in Section 
3.4, a single jump is not enough). The Malliavin calculus conveniently enables us to bypass intractable 
direct estimate of the characteristic function of L(T~ 1 I' 1 Ht), and results in fairly simple conditions. 
Our result has the following statistical implication. Suppose that we can directly observe {X t : < 
t < T}, based on which we want to estimate 6>o := -E[Ao] (the mean of the stationary distribution). 
A natural estimator is then given by 

1 ' T 



9 T := - / X s ds 
1 Jo 

We easily see that T~ x / 2 Ht = T x / 2 {6t — 6q) with j3 — 1 and 7 = p = 0, hence the consistency, 
asymptotic normality, and higher order expansion of 6t are obtained according to our result. 

The main result is given in Section 2, followed by the proof in Section 3. 

2 Edgeworth expansion 

2.1 Statement of result 

We are given a stochastic basis (Q, J 7 , F = (Ft)teM+,P), on which our processes are defined. 

Assumption 2.1. X is strictly stationary with a stationary distribution F <= |~) >Q L P (P). 

We remark that: under Assumption 2.1 X is exponentially /3-mixing and ergodic; Assumption 2.1 
is equivalent to Z\ € Pl P >o L P {P)- See [4] for details. 

Denote by (6, C, II) the generating triplet of Z in the form 

ip(u;Z t ) =explt(ibu- \-Cu 2 + / (e mz - 1 - iuz)Ii(dz) J i, 

where b £ R, C > 0, and the Levy measure II defined on R is a cr-finite measure satisfying II({0}) = 
and Jn<l Ki z2 n(<iz) < 00. Then the process H of (3) satisfies 

dH t =P(X t - 4 1} )dt + pdZ t , H = 0, 

where Z t := Z t - E[Z t ] = Z t - E[Z x ]t and 

nf ] :=i -fe d„log-E[exp(m£)], 

the fc-th cumulant of £, with d v denoting the (partial) differentiation with respect to a variable v. 
Denote by A the Poisson random measure associated with jumps of Z. We decompose it as 

A(dt, dz) = p (dt, dz) + p(dt, dz) 

for some Poisson random measures fj, and /1; by the independently scattered property of A, such a 
decomposition is always possible. Correspondingly, we write 

tt(dz) = v h {dz) + v(dz), 

where v and v stand for the Levy measures on R + associated with /r and p, respectively. 

Assumption 2.2. Either one of the following two conditions holds true: 

(i) C > (no condition is imposed on the jump-part characteristic) ; 



(ii) C — and there exists a non-empty open subset of M\{0} on which v admits a positive C 3 - 
density, say g, with respect to the Lebesgue measure. 

Let us introduce the notation necessary for the Edgeworth expansion; see [6] for more details. We 
introduce the r-th cumulant function of T~ x I 2 Ht (r 6 N, r > 2): 

Xr,r(«) := d r u logS [exp^T- 1 / 2 ^) 

Let p > 3 be an integer. The {p — 2)-th Edgeworth expansion ^ Pj t (a signed measure) is defined by 
the Fourier inversion of u n- ^ Pi t(m), where 



/i \ p- 2 

$ p , T (u) := exp - X T,a(«) + ^T~'"/ 2 P T , T (iO, 



r=l 



with P t ^t{u) specified via the formal expansion 



r\ 

r=2 



Let (/)(•; E) stand for the one-dimensional centered normal density having variance E > 0, then the 
r-th Hcrmite polynomial associated with (/>(■; E) is h r (y; E) := (— l) r cf)(y; E) _1 dy<f>(y; E). Let 



Xr.T ■= (-l) r Xr,T(0), 

the r-th cumulant of T-^ 2 H T ; in Section 2.2, we will see that Xt.t = O^-^-^/ 2 ) as T -> oo. The 
density of ^ p ,t with respect to the Lebesgue measure is given by 



ffp 



(y;T-V 2 H T ) 



fe=l j=i «n fc,eN: v i / v i i ; 

felH hfc,=fc 



where E T := X 2,t; wc will approximate ElffT-^Hr)] by %, T [f] := / f(y)g P (y;T- 1 / 2 H T )dy. 

Let po : = 2[p/2] and denote by £(M,po) the set of all measurable functions / : R — > R satisfying 
\f(x)\ < M(l + |a;|P ) for every a; £ R. 

Now we can state the main result. 

Theorem 2.3. Let X,Y,H be given through (1) and (3), and suppose that (2) and Assumptions 2.1 
and 2.2 hold true. Fix any positive number E° such that 

E°>^(/3 + pA) 2 4 2) . 

Then, for any M, K > 0, there exist positive constants M* and S* such that 

E[f(T-^ 2 H T )} - %, T [f] < M* f sup \f(x + y)- f(x)\<f>(x; E°)dx + (T-&- a +*'>/ a ) (4) 

JR \y\<T-'< 
for T — > oo uniformly in f £ £{M,pa). 

Most often in practice, the first term in the upper bound in (4) can be quickly vanishing by taking 
K large; for example, it is the case when / is an indicator function f = 1a for various icK, such as 

A = (— oo, a], A — [a, b], and so on. 

2.2 Explicit coefficients 

The approximating density g p (-; T~ 1 > 2 Ht) involves the cumulants X2,t, X3,t, ■ ■ ■ , Xp,t- We here prove 
the explicit formula for Xr,T, r > 2. 



Noticing the explicit solution X t = e xt Xo + J Q e A '* s ^dZ s , we can apply the stochastic Fubini 
theorem to obtain the relation 



X s ds = 77 (A, t)X + / ?/(A, t - s)dZ s 



(5) 



where r](X, u) = A x (l — e Au ); one can consults [2] for a detailed analysis of integrated OU processes, 
especially in the context of financial econometrics. It follows from (1), (5), and the special relation 

kXn F = K z for fc G N (see [1, 4]) that we can express Ht as 

H T =f3 V (\,T)X Q ~T(P + p\)4 ) + ( { P + f3 V (\,T~ S )}dZ s . 



Hence, using the independence between Xq and Z we obtain 

r'-l' 



Xr.T = {-if 



()' 



Ik (pr-WriiX, T)u; f) + j 5> ({p + ^(A, T - s)}T- l ' 2 u- Z x ) ds 



u=0 



{/3T- 1 /2 r?(A) T)} r 4 r) + f T (ip + v (X,v)}T-'/ 2 y 'dvXrnP 



ji-(r-2)/2 



T- 1 {pr,(X, T)} r + XrT- 1 / {p + f3 V (X, v)} r dv 

Jo 



„ F . 



By making use of the differential equation 9 S {?y(A, s)} = k {rj(X, s)} — Xk {r)(X, s)} with rj(X, 0) = 
and then integrating the both sides with respect to s over [0, T], we can proceed as in [5, Section 3] 
to conclude that 



rr—(r —21/2 
Xr.T = T y I' 

where M. t ,t{J) is given by 

■M r ,r(0) = 1, 



T- 1 {pr,(X, T)} r + Xr Y, I ) P r - j P j M r , T (j) 

3=0 V/ 






(6) 



M-.rtt) - A-J - T- 1 A-(^+ 1 ) J2 k ^ iM\ T)} k , 3 > 1. 



fc=i 



Thus we can explicitly write down the coefficients of the Edgcworth expansion ^ p _t up to any order. 
It is obvious from (6) that X t.t = OiT"^- 2 ^ 2 ) for r > 2; 



TC- 2 )/V,T^Ar^ r W^A 



j=o 



-J K M 



In particular, 



Sr = X2.r^2A- 1 (/3 + pA) 2 4 2) , 
hence the necessity of the condition (2). 



3 Proof of Theorem 2.3 

We will apply [6, Theorem 1]. In order to ensure distributional regularity necessary for the Edgeworth 
expansion, we will make use of a Malliavin calculus with an effective truncation functional. The main 
idea of the proof is in principle similar to that of [5, Section 4] treating the stochastic volatility model, 
where X expresses the latent positive volatility process. However, the OU process X in the present 
model can take negative values too, and, as such, the way of constructing a truncation functional 
is essentially different from that of [5]. To save space, we will sometimes omit the technical details, 
referring to the pertinent parts of [3, 5]. 

Let us briefly overview the fundamental device. By means of [6, Theorem 1], it suffices to verify 
the following conditions: 



[Al] X is strongly mixing with exponential rate; 
L42] sup te[0T ] ||ift|| L p + i( P ) < oo for each T <G R+; 

[A3] there exist positive constants t , a, a 1 and B, and a truncation functional ip : (Q,J-) 
([0, 1], B([0, 1])) such that < a, a' < 1, 4a' < (a - l) 2 anc 



< a , 



E sup \E[ipex.p(iuH t o)\X ,X t o] 

. \u\>B 

1 - E[tp] < a. 

As was mentioned in Section 2, Assumption 2.1 ensures [Al] and [A2] (see (6)), so that it remains 
to verify [A3], which is a version of conditional Cramer conditions. Although it may be difficult in 
general to verify [A3] , we will be able to construct a specific truncation ip which significantly simplify 
the task. 

We also note that the condition (A' — 4) of [3, p. 60 and p. 130] (smoothness of the coefficients, 
and intcgrability under cut-off through an auxiliary function) is indispensable. We will mention this 
point in Section 3.2 

3.1 Transformation of the Poisson random measure 

In order to execute a Malliavin calculus of [3] , we introduce a transformation of the absolutely con- 
tinuous part of the Poisson random measure. 

Under Assumption 2.1, Z admits a Levy-Ito decomposition of the form 

Z t = \n { F ] t + VCw t + zfj,\ds, dz)+ zjx(ds, dz), t e K + , 

JoJs. JoJr 

where w stands for a one-dimensional Wiener process defined on (fi,.F, F,P), fi b (dt,dz) := fi b (dt,dz) — 
v b (dz)dt, and fl(dt,dz) := fj,(dt,dz) — v{dz)dt. 

Assumption 2.2 assures the existence of a bounded domain 

£o = (ci,c 2 )cl\{0}, 

for which the Levy density g of v satisfies that 

inf g(z) > 0. 

z£E 

Without loss of generality, we may and do suppose that Ci,C2 > 0: if ^(K + ) = 0, then take -Z as 
Z anew. We introduce the change of variables z* — z*{z) = g + {z) through z* = z*{z) — J g(v)dv 
for z € £o; obviously, g + is strictly decreasing on Eq. Let g~ denote the strictly decreasing inverse 
function of g + defined on 

£ = (.g+( C2 ),.g+( Cl )). 

Let n* denote the integer-valued random measure defined by 

t pa2 pt pg (ai) 

h(s , z) /i(ds , dz) = / / h(s,g~(z*))/j,*(ds,dz*) 

0Ja t JoJg+(a 2 ) 

for each t S R+ , 01,02 € R such that a\ < ai , and for any measurable function h on R+ x R + ; in 
particular, 

E\p*{[0,t],B)] =thch(B). 

Writing fj,*(dt,dz*) = /j,*(dt,dz*) — dtdz*, we transform fj, (on [0,i] x Eq) into fj,* as follows: 

rt pc 2 pt pg + (c 1 ) 

/ / zfx(ds,dz)= / / g-{z*)fi*{ds,dz*). 

JoJct JoJg+(c 2 ) 



The bivariate process (X, H) satisfies the stochastic differential equation 

;s)=<^.>(>+< 

z( 1 )((j b + l Es fi)(dt,dz)+ f J(z*)( 1 )(i*(dt,dz*), 

\Pj JEL>[g+( Cl ),ac) \PJ 

where J{z*) := g~{z*)l E {z*) for z* e PU[g + (ci),oo). As g~ is strictly decreasing, we have \dJ(z*)\ > 
for z* e EU [g + (d),oo). 

3.2 Malliavin covariance matrix 

Fix any constant t° > and define (Cl, B, P) to be the Wiener-Poisson canonical space (see [5, the last 
paragraph in page 1178]), on which we are given the flow (X (• , v) , H (• , v)) T associated with (X,H) 
of (7) starting from v = (x, h) T € R 2 : 

X(MA - fA + [\n^-X(s,v))^y s + VcQw t 



H(t, v )j \hj y * "\-PJ vv 

t r /1\ ft 

z( 

Pj Jo ^BU[g+(ci),oo) 



1 )(^ + l Es fi)(ds,dz)+ [ [ J(z*)r)ji*{ds,dz*). 

PJ h JEUfo+( Cl ),oo) \PJ 



Under the present assumption, the flow (X(-, v), H(-, v)) T clearly satisfies the condition (A' — 4). 

Let a; be a random variable independent of (w, fi +l E cfj,, (i*) such that C(x\P) = F (the distribution 
under P), and v := (x, 0) T . We will compute the Malliavin covariance matrix of (X(t°, v), H(t°, v)) T , 
whose "non-degeneracy" is essential here. 

Let QeI 2 ®K 2 bc given by 



P y 

In view of (7), the process K(t,v) := d v (X(t,v),H(t,v)) T satisfies that, for each v, 

d rs„ ^ ( -\d x X(t,v) \ ~ T , I± , 



so that 

(Different from [5, Eq.(29) in page 1181], S(-,v) is independent of v.) 

Pick positive constants c'- and c" (j = 1,2) in such a way that < c\ < c[ < c'{ < c'J, < c' 2 < c 2 < 
oo, and let 

E:=(g+(4),g+(c'{)). 

(Trivially, E <s E.) Let r\ : R + — > E + be any bounded smooth function satisfying the conditions: 

(i) ' m i z . eE rj(z*) > 0; 

(ii) r,(z*) = 0{oiz*t(g+(c' 2 ),g+(c' 1 )). 
The Malliavin covariance matrix of (X(t°, v), H(t°, v)) T is then well-defined and given by 

U(t°,v) := exp{t°Q)S{t°,v)exp{t°Q T ), 
where 

S(t,v) = C J cxp(-sQ) f 1 P A cxp(-sQ T )ds 

+ / / exp(-sQ) ( l P 2 ) cxp(-sQ T )l/(z*)/i*(ds,dz*), (8) 

JoJ-E \P P J 

with V(z*) := {3J(z*)} 2 ?7(z*); see [3, Section 10] for details of (8). Thus we arrive at the identity 

detU(t ,v)=e~ 2Xt "detS{t°,v), a.s. (9) 



3.3 Completion of the proof under Assumption 2.2 (i) 

Suppose that C > 0. It follows from (8) that, in the matrix sense, 

ft 



+0 „-.\ ^ r< I „-sQ I 1 P 



H2 sym. 

X Hi - W/X)H 2 x 2 t° - 2(p/X) X H 1 + (P/XfH, 

,0 

where Hk := J Q e kXs ds and \ :== P + /^/A. The determinant of the rightmost side is 

C 2 \- 4 ((3 + p\) 2 l^(e 2Xt " - 1) - (e At ° - l) s 

which is positive as soon as t°X 7^ and j3 + pX 7^ 0. Thus S(t ,v) is bounded from below by a 
positive-definite matrix, hence the non-degeneracy of U(t°, v) follows from (9) without any non-trivial 
truncation functional; simply let i/> = 1 in [A3]. Thus we have obtained the non-degeneracy of the 
Malliavin covariance matrix (i.e. enough integrability of {detU(t°, w)} _1 ), which corresponds to [5, 
Lemma 6]. 

We further notice the following. 

• The flow (X(t,v),H(t, v))te[o,t°] satisfies the condition (A' — 4) (as was seen in Section 3.2), 
hence the analogous assertions as [5, Lemmas 7] holds true. 

• Following the same argument as in [5, pp. 1184-1185], we see that there exists a random variable 
$£0 G L X {P) such that 



E 



sup \E[cxp(iuH t o)\Xo,X t o] 

\u\>B 



< jjE[\&* 



for every B > 0. 



After all, we have deduced the analogous assertions to [5, Lemmas 6, 7 and 8], completing the 
proof of Theorem 2.3 under Assumption 2.1 and Assumption 2.2 (i). 

3.4 Construction of a truncation functional 

It remains to prove Theorem 2.3 under Assumptions 2.1 and 2.2 (ii); then, in order to verify dis- 
tributional regularity we have to make an effective use of jumps. We will construct the truncation 
functional ip in an explicit way through two diffusive jumps. 

We continue the argument of Section 3.2. Let £1,^2 G (0,i°) be constants such that t\ < t-2, and 
fix Zq G E. Let e > be sufficiently small so that: 

. If nTf = for p. := (tj - e,tj + e), 3 = 1,2; 
• 9 + (4) <z -e<z Q + e< g+(c'{). 
Let E e := (zq — e, zq + e) and 

A*:={»*(I*,E e ) = l fori = 1,2.}. (10) 

According to the independently scattered property of p* and since the Levy measure associated with 
p* over E) here is the Lebesgue one, we have 

P[A f -} = {P[p* ([0,26], [0,2e])]} 2 = {4 e 2 exp(-4e 2 )} 2 > 

for each e > 0. 

We define the truncation functional ijj e by ^ = C(^e); where ( : M + — > [0, 1] is a non-increasing 
smooth function such that £(x) = 1 if < x < 1/2 and ((x) = if x > 1, where 

2 
^ e = 1 + 3dctC/(iO, y) ■ (11) 



We will show that the Malliavin covariance matrix U(t , v) is non-degenerate on the event A e for any 
e > small enough. 
Noting that 



sup 

s:\s-tj | <€ 



-sQ 



pX-^l-e^) 1 



+ sup \V{z) - V(z )\ -> 

z: \z — zq I <e 



as e — > and by virtue of (10), we apply Taylor's expansion around zq and tj (j = 1,2) on A e to 
conclude that 



s( t °,v)>Y: 



j=l «I*JE< 



1 p 

p P 2 



V{z*)fi*{ds,dz*) 



j~i \P P 

V(z )M' + o(l) 



e-^ Q V(z ) + o(l) 



as e — ¥ (we use the symbol o(l) for matrices too), where 



M e 



with J( ] 



J(2) 



sym. 



(p + (iX- 1 ) J« - /3A" 1 JW 20 + /3A- 1 ) 2 - 2/3A- 1 (p + /SA^ 1 ) i« + /3 2 A~ 2 J^ 



D \t ± 



e A ' 2 and j( 2 ) := e 2Atl + e 2A * 2 . Therefore 

detS(i°, v) > V{z a f\- 2 {P + Ap) 2 (e Atl - e A ' 2 ) 2 + o(l), 

which is positive for e sufficiently small whenever pX + (3 =/= and ti =/= t 2 - [We note that a single jump 
is not enough: if we instead estimate S(t°,v) as 



°,*)>e-^Q £) e -^V(z ) + O (l), 



S(t 



then the determinant of the first term in the right-hand side turns out to be identically 0.] 

We may set V{zq) arbitrarily large by choosing the function r\ suitably, so that, recalling (9) we 
conclude that detZ7(£° ,v) > 1 on A e for some e > 0. Fix such such 77 and e > 0. The definition (11) 
leads to the estimate 

P[i < 1/2] > P [{detU(t°,v) > 1} n A e ] = P[A € ] > 0, 

hence the assertion corresponding to [5, Lemma 6] holds true. We fix the r\ and e > in the sequel. 
Clearly i/j e > implies that 1/3 < detU(t°,v), hence 

^{dctuit ,^ 1 e f) LP(P). 

0<p<oo 

This implies that the integration-by-parts formula under the truncation ip e is in force. Then, as before, 
we could deduce the assertions corresponding to [5, Lemmas 7 and 8]: 

• The flow (X(t, v), H(t,v)) te [ 0t o] satisfies the condition (A' — 4) (as was seen in Section 3.2); 

• There exists a random variable $" G L l (P) such that 



E 



sup 

|u|>B 



E[ip e exp(iuH t o)\X ,X t o] 



1 
B' 



't°\ 



for every B > 0. 
The proof of Theorem 2.3 is thus complete. 



Acknowledgement . 

This work was partly supported by JSPS KAKENHI Grant Number 23740082 (H. Masuda). 

References 

[1] Barndorff- Nielsen, O. E. and Shephard, N. (2001), Non-Gaussian OU based models and some of 
their uses in financial economics (with discussion). J. R. Stat. Soc. Ser. B Stat. Methodol. 63, 
167-241. 

[2] Barndorff-Nielsen, O. E. and Shephard, N. (2003), Integrated OU processes and non-Gaussian 
OU-based stochastic volatility models, Scand. J. of Statist. 30, 277-295. 

[3] Bichtclcr, K., Gravcreaux, J. and Jacod, J. (1987), Malliavin calculus for processes with jumps. 
Stochastics Monographs, 2. Gordon and Breach Science Publishers, New York. 

[4] Masuda, H. (2004), On multidimensional Ornstein-Uhlenbeck processes driven by a general Levy 
process. Bernoulli 10, 97-120. 

[5] Masuda, H. and Yoshida, N. (2005), Asymptotic expansion for Barndorff-Nielsen and Shcphard's 
stochastic volatility model. Stochastic Process. Appl. 115, 1167-1186. 

[6] Yoshida, N. (2004), Partial mixing and Edgeworth expansion. Probab. Theory Related Fields 129, 
559-624. 



